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Abstract 

Usually, neutrino oscillation experiments are analyzed within the two-flavor frame- 
work which is governed by 1 mass-squared difference and 1 mixing angle. But there 
are 6 parameters, 2 mass-squared differences, 3 mixing angles, and 1 CP phase within 
the three-flavor framework. In this article, we estimate the effect from the smaller 
mass-squared difference, the other mixing angles, and the CP phase, which we call 
three- flavor effect, for the determination of the mass-squared difference and the mix- 
ing angle from the z/^'s survival and transition probability with the two- flavor analy- 
sis. It is found that the mass-squared difference from the two-flavor analysis is slightly 
shifted from the larger mass-squared difference by the three- flavor effect. The order 
of magnitude of the three-flavor effect for the mass-squared difference is compara- 
ble with that of the expected error for the mass-squared difference of the two-flavor 
analysis in the future long base-line neutrino oscillation experiments. The CP phase 
dependence of the — > Ue transition probability is also shown. 



1 introduction 

The three neutrino framework has 9 physical parameters: 3 neutrino masses, 3 mixing 
angles, and 3 CP violating phase, if neutrinos are Majorana particles. Neutrino oscillation 
experiments are sensitive to 6 parameters: 2 mass-squared difference, 3 mixing angles, and 
1 CP phase. Usually, data from the experiments are analyzed within the two-flavor frame- 
work, which is governed by only 1 mass-squared difference and 1 mixing angle. So far, for 
the long base-line neutrino oscillation experiments, we have been able to neglect the effect 
from the smaller mass-squared difference, the other mixing angles, and CP phase, which 
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we call three-flavor effect. This is because the error of the larger mass-squared difference, 
related to the atmospheric neutrino observations [1] and the on-going long base-line exper- 
iment K2K [2], is larger than the smaller mass-squared difference, which is obtained from 
solar neutrino observations [3] and the reactor neutrino experiment KamLAND [4]. 

The long base-line neutrino oscillation experiments in future [5, 6, 7] plan to measure the 
mass-squared difference and mixing parameter precisely. Because the order of magnitude of 
the ratio between the smaller mass-squared difference and the larger one is supposed to be 
similar to that of expected error of the future long base-line experiments, it is necessary to 
take into account the contribution of the three-flavor effect in the determination of the mass- 
squared difference and mixing angle from the future long base-line precision measurements. 
In this article, we estimate this effect using the three-ffavor framework. We discard the 
result of the LSND experiment [8]. Obviously, this analysis is, in general, not valid, if the 
LSND result is confirmed by the MiniBooNE experiment [9]. From the survival probability 
of i/^, we find that the larger mass-squared difference is shifted by the three-fiavor effect and 
that the order of magnitude of this shift depends on the neutrino energy and is similar to 
that of the smaller mass-squared difference. We obtain the same result from the transition 
probability i'^ ^ phase dependence for the transition probability. 

A lot of groups have analyzed the experimental data with the three-flavor framework 
numerically [10]. However, these analyses cannot point out the speciflc reason for the 
value of parameters. In this paper, we point out the speciflc contribution to the parameters 
from the three-flavor effect. We think that these formulations are useful to study in the 
numerical analysis. 

This article is organized as follows. In section 2, we will show the useful notations and 
the convenient form of the probability for easy estimating the contribution of the three- 
flavor effect. In section 3, we will estimate the three-flavor effect for the i/^ disappearance 
mode. We will also estimate the three-flavor effect for the v^j^ — > transition mode, in 
section 4. Finally, we will be devoted to the summary in the last section. 

2 notations 

In the three neutrino framework, and in the basis in which the charged leptons are diagonal, 
the three weak interaction eigenstates, {a — e, r) are expressed as 

3 

= I](^MNs)m , (1) 
i=l 

where are the three mass eigenstates and Vmns is the Maki-Nakagawa-Sakata (MNS) 
matrix [11]. We adopt the following parameterization [12] 

Vrnm = UV = U di^g{e^'^^'\ e'^'^'\ 1) , (2) 



where V cannot be determined from the neutrino oscillation experiment. The matrix U, 
which has three mixing angles and one phase, can be parameterized in the same way as 
the Cabibbo-Kobayashi-Maskawa matrix [13]. Because the present neutrino oscillation 
experiments constrain directly the elements Uf,2) Ue3, and we find it most convenient 
to adopt the parameterization [14, 15], where these three elements in the upper-right corner 
of the U matrix are the independent parameters. Without losing generality, we can take 
Ue2 and to be real and non-negative. By allowing 1/^3 to have the complex phase, 

C/e2,C/^3>0, Ues=\Ue3\e-''^ (0 < < 27r) , (3) 

these Ue2, Uij,^, [C/esI, and (f) are the four independent parameters. All the other elements of 
U are then determined by the unitary conditions, 



U,, = Jl- |t/e3|2 -\U,2\\ U,s = Jl- |t/e3|2 - \U^-,\^ , (4a) 



Ue2U,3 UMlUk jr _ U^^Uel + ^e2^.3^e*3 ^ 



For the convenience, the independent parameters in the MNS matrix are rewritten as 

C/g3 = sin 6*13 , C/g2 = sin 6*12 cos 6*13 , C/^3 = sin 6*23 cos 6*13 . (5) 

The atmospheric neutrino oscillation experiments, which measure the survival prob- 
ability determine the absolute values of the larger mass-squared differences and one-mixing 
angle [1] as 

< 3.4 ± 0.5 X lO-^^eV^ , and sin^ 2^atm > 0.92 , (6) 

at the 90% C.L. The K2K experiment [2] confirms the above results. These values are 
planed to measure more precisely, a few percent order by the future long base-line experi- 
ments [5, 6, 7]. The solar neutrino experiments, which measure the Ve survival probability 
in the sun [3], and the KamLAND experiment which measure the survival probability 
from the reactors [4], determine the smaller mass-squared difference and another mixing 
angle as 

(5rr4i = 8.2+°:^ x 10"^eV=^ , and tan^ ^^oi = 0.40+[!:J^? . (7) 

It is remarkable point that the order of the smaller mass squared-difference is as same as 
that of expected error of the future long base-line experiments [5, 6, 7]. Thus, it is necessary 
to take into account the contribution of the three-flavor effect in the determination of 
the mass-squared difference and mixing angle analytically, because experimental data are 
analyzed ordinary in the two-flavor framework. The CHOOZ reactor experiment [16] gives 



1.5 X lO-'^ < 



the upper bound of the third mixing angle as 



sin^ 2^,ct < 0.20 for 6m^ = 2.0 x lO^'^eV^ , 
sin^ 2^,ct < 0.16 for W = 2.5 x lO^-'^cV^ , 
sin^ 29,a < 0.14 for dm^ = 3.0 x lO'^eV^ , 



(8) 



at the 90% C.L. Since we can always take |Ai2| < |Ai3| without loosing generality, we 
assume that A 12 is from the results of solar neutrino and reactor anti-neutrino observations 
and Ai3 is related to the atmospheric and long base-line neutrino experiments. The sign 
of A12 is determined form the solar neutrino observation, A12 > 0. However, that of 
Ai3 cannot be determined by any observations. In this article, we call A13 > "normal 
hierarchy" and A13 < "inverted hierarchy". Under the |Ai2| < |Ai3| relation, ^atm, ^soi, 
and 9^ct a-re related to the MNS matrix elements as : 



2 It/, 



e3l 



1 - \/l -sin^2^. 



ret ) 



i-i^esr- 



:i-|t^e3P) -sin^2^soi- 



(9) 



Starting from the flavor eigenstate a, the probability of finding the flavor eigenstate /3 
at the base-line length L is, in vacuum. 



3 / ^2 \ 

E(^MNs)/3jexp {VuNsTaj 



(10a) 



6^^ - 4Re {u^sU;,U^,U:, + U^,U;,U0,U:,} sin 



2 ^13 



+ 2Re {u^2U;2U^3Ks) (sin A12 sin A13 + 4 sin^ ^ sin^ ^ 



4jK(-n^^sinA,2-sin^^sinA 



13 , 



(10b) 



where Jum the Jarlskog parameter [17]: 



J. 



MNS 



Im ( ( ^MNS ) ai ( ^MNS ) ( ^MNS ) ft' ( VuNS ) *aj ) 

-I- — K^f!3 



= A sin , 



(11) 



which is defined to be positive for {a, /?) = (e, //), {/i, r), (r, e) and — (1, 2), (2, 3), (3, 1). 
A is the absolute value of the Jarlskog parameter. In addition, A^j is 



2 2 

— mf 



Sm 



2E„ 



2E, 



'^Lc:^ 2.534- 



L(km) , 



(12) 



where Ei, is the neutrino energy. 

We rewrite eq. (10) in a form convenient to estimate the contribution of the smaller 
mass-squared difference, 



P,^_„ = Po(a, P) + Pi (a, P) X sin A^^ + P^ia, /3) x 4 sin 



2^12 

2 ' 



where each term is 



2 ^13 



2 ' 



Po(a, (5) = Sap - 4Re {u^sU^.U^.K, + t/„2^^2^/33^;3} sin^ 

Pi(a, (3) = 2Re (t/«2t^;2t^/33t^:3) sin A,, - 4 J^^ sin^ ^ , 

P2(a, (3) = -Re {u,,U;,Up,U:, + U^^U^p^U^sK, cos A13) + jiSfg^ sin A 



13 



(13) 

(14a) 

(14b) 
(14c) 



All the above formulas remain valid when replacing the mass-squared differences and 
the MNS matrix elements with the "effective" ones. 



(15) 



as long as the matter density remains the same along the base-hne. The effective parameters 
Uai are defined from the following Hamiltonian 
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(16) 



and Ajj is defined as 



A,-, = —^L = l„ L . 



2E^ 2E^ 
The term a in eq. (16) stands for the matter effect [18], 

a(E^) = 2V2GFneE^ = 7.56 x 10-^(eV2) ^ 



g/cmV vGeV 



E,. 



(17) 



(18) 



where ng is the electron density of the matter, Gp is the Fermi constant, and p is the matter 
density. When Smf2 < a < Sml^ and C/^g <^ 0(1), the effective mass-squared differences 
are written as 



V13 ~ - cos , ^12 - - cos 2^12 , 



and the mixing angles become 



(19) 



^23-^23, 



13 



^1^,3, tan 2^, 



(5m^2 sin 2^;^2 
(5mf2COs2^i9 — a 



12 



(20) 



Hereafter, we use U^^ and A^j instead of C/^^ and Aj^, because of simplicity. But we have 
to keep in mind that these values depend on the neutrino energy. 



3 survival probability of u 



Prom eqs. (13) and (14), the survival probabihty of is written as 



Pu^^u^ = Po{l^,l^) + A (a*, A*) X sin Ai2 + P2{p, 1^) x 4sin^ ^ , 



where 



Po(a^,A^) = 1-4 



Pi(//,//) = 2 



P2{f^,fJ') = 



u. 



H2 



/i3 



u. 



u. 



/u3 



2\ . 2 ^13 



sm 



2 ' 



U. 



H2 



/u3 



sin Ai3 , 



Ml 



+ 



/i3 



COS A 



13 



The survival probability of u in the two-flavor framework is written as 



(21) 

(22a) 
(22b) 
(22c) 



p(2) 



1 -sin2 2e„„sin2^. 



(23) 



where 6^^ is the mixing angle and A^^ is the mass-squared difference, which are obtained 
from the two-flavor analysis. We expect that all these parameters, especially A^^, to be 
shifted by the three-flavor effect. In order to estimate this effect, we rewrite them as 



sin^ 29^ 



1.0 -£, 



Ai3 + 25, 



(24) 



where 5 denotes the three-flavor effect, and e stands for the difference from the maximal 
mixing. Both of them depend on the neutrino energy, base-line length, and the oscillation 
parameters. We assume that 5 is smaller than A13 and A13 ~ 0(1) for long base-line 
experiments. 



Prom the atmospheric neutrino observations and K2K experiment, wc aheady know 
that oscillate to another flavor maximally at the first-dip, A^^{E^ — £^dip) — for 
normal hierarchy. By using eq. (24), i.e., A13 — n — 2(5^^^, eq. (22) becomes 



Po(/i,/i) = l-4 



cos' 



Pi(//,/i) = 2 



US'' Ut^'sm25''^, 



u: 



dip 



u: 



dip 
111 



dip 
fj,3 



COS 



and eq. (23) 



p(2) 



^dip 



(25a) 
(25b) 
(25c) 

(26) 



where the label "dip" in \U^^\, d, and e means that these quantities take some fixed value 
at the first-dip energy £^dip- From eqs. (25) and (26), we obtain 



^dip 



Ut^\"{l-\ui^\")-{l-e'^^) 1 



u: 



dip 



/i2 



dip 
/i3 



u: 



dip 



/i3 



(27) 



The first term of eq. (27) has to be zero because of the assumption 0{d) < 1, and therefore 
we obtain 



2 (1 ± 



{21 



When we take a negative sign in eq. (28), t^'^'P becomes 



A dip 



ry-dip 



+ 1- 



rrdip 



(29) 



Since the best fit value of the mixing angle is maximum from the experiments [1, 2], we 
take £^'P = 0. Thus, simphfies to 



A^^P 



ry-dip 



Prom eq. (24), the larger mass-squared difference at E'dip is 



where 



2 _ 27ri?dip 
<^^dip = — > 



(30) 

(31) 

(32) 



is from the two-flavor analysis. Equation (31) denotes that the order of magnitude of the 
three-flavor effect is roughly the same as that of the expected error of future experiments. 



Since we know that both 5m^2 ^^^1 



dip 



ii2 



are positive, the mass-squared difference from 
two-flavor analysis is slightly smaller than the larger mass-squared difference. Wc also 
obtain the relation between 5 and the smaller mass-squared difference at the first peak 
-£^peak, where A^^(£;peak) = 27r, 



^peak ^ 



^peak 



1 - 



^peak 



(33) 



Thus, the larger mass-squared difference at E^^ak is 



(5mi3(£;peak) = 5ml^^ + 5m?2(-Epeak) 



peak 



jJL2 



where 



<^"^peak 



47r£'peak 



(34) 



(35) 



is also from the two-flavor analysis. We obtain the same results for the inverted hierarchy. 
Since 5m{^ is not changed by matter effect at Ei, <10 GeV, we obtain the relation between 
eq. (31) and eq. (34): 



5m^ip + (5m?2(-Bdip) 



^.dip 



^"^peak + (^mi2(£^peaic) 1 - 



r 7-pcak 



(36) 



By using the deflnition of 5^^dip,peak' '^^ 



L 



j-jdip 



T /-peak 
,,2 



2n 



2E, 



peak 



E, 



dip 



±1, 



(37) 



where the sign of the right-hand side corresponds to the type of the mass hierarchy, the 
positive sign being that of the normal hierarchy. Prom the eqs. (31) and (34), we easily 
understand that the three-flavor effect for the larger mass-squared difference depends on 
the energy. When p = 2.5(g/cm^) and E^^ ~ 0(l)GeV, is to shift away from 6^2 towards 
90°, which is obtained from eq. (19). From eqs. (4) and (5), the value of \U^2\ becomes 0. 
Because Sml2{Eu) is also changed by the matter effect, which is shown in eq. (19), eq. (37) 
becomes 



L a{Epeak) 



8m\2 cos 29-^2 



27r 



-E^dip — 2E, 



±1 



peaJf 



(3^ 



where the value of 5m\2 and 6^2 is that of vacuum one, which are listed in eq. (7). This 
relation suggest us that we can determine the sign of the 5m\'^ from the long base-line 
experiments, when we measure the energy of "peak" and "dip" precisely and we know the 



value of the smaller mass-squared difference, mixing angle with small errors. This result 
cannot be obtained from the numerical analysis. This fact points out that we can pick up 
the three-flavor effect from the fltting function of the survival probability which is obtained 
from the experimental data. 

4 transition probability 

From eqs. (10) and (14), the transition probability, — > is written as 



where 



Pv^^v^ = Poil^, e) + Pi(/i, e) X sin A12 + P2{n, e) x 4 sin^ ^ , 



-Po(/^,e) 



sm 



2 Al3 

2 ' 



Pi(//, e) = 2 (Re {u:,U^,U,,U;,) sin A13 + 2 J&=^ sin 



2 ^13 



P2(//, e) = -Re (U^,U:,U,,U;, + U^^Ul^U^^U;^ COS A13) - J&1 sin A13 ■ 



(39) 

(40a) 

(40b) 
(40c) 



Here, Jmns — — ^sin^ = — J. Under the two-flavor framework, this transition probability, 
— > i/g is written as 



^i:U=sin2^,.sin^^, (41) 

where A^g is the mass-squared difference and is the unknown mixing angle. We suppose 
that these two parameter, especially A^e, are changed by the three-flavor effect. As in the 
case of the survival probability, we rewrite these parameters as 



sin^ 6l„e = ^ , 



A13 + 25, 



(42) 



where S is smaller than A13 and h, in general, can take an arbitrary value. Before estimating 
the three-flavor effect, let us calculate the value of h and S for A12 = 0. The transition 
probability becomes 



P —A 



2 . 2^13 

sm 



2 ' 



(43) 



with A12 = 0. At the flrst peak of transition probability, A^e(-E'peak) = tt, eqs. (41) and (43) 
become 



P„(A^e = 7r) = 4 



cos^ S , 



(44) 



respectively. Prom these equations, h and S are solved as 



= /i(Ai2 = 0, A^e = tt) = 4 
5o = 5(Ai2 = 0,A^e = vr) = 0, 



when we suppose that A^g is same as A13 without three-flavor effect. 
When we assume that A12 is nonvanishing, but small, eq. (39) becomes 

P,^^,^ = Po(//, e) + Pi(Ai, e) X A12 + P2(//, e) x Aj^ + 0{Al^) , 

where Po(/x, e), Pi{iJ,,e) and P2(/x, e) are 



Po(//,e)=4 
= 4 



/i3 



|C/g3|^ COS^ 5 



/i3 



Pi(//, e) = 2Re^23 ^in 25 - 4 Jcos^ (5 
= 45Re^23-4J + 0(52), 

P2(/i, e) = — Re^]^2 + I^c^23 25 + J sin 26 
= -Re>li2 + Re^23 + 25 J + 0(5^) , 

at the first peak Ep^ak- Here, the symbols A,^, are defined as 



A.. — TJ .TJ*.TJ .TJ* . 



(45) 



(46) 



(47a) 

(47b) 
(47c) 

(48) 



The value of U^j, A^j, and 5 is fixed at ii^peak*- We assume that h at ii^peak is function of 
A12, 



where ho is equal to 4 U^^U^^ 
eqs. (47) and (49), 5 can be so 



and are independent of the neutrino energy. Prom 
ved as 



(ai + 4 J) + (Re^i2 " Re^23 + «2) A12 + a3Af2 



4Re>l23 + 2 JA 



12 



(50) 



When we take the limit A12 0, S must vanish because of eq. (45). Thus, 



and 



= -4J, 



^ _ (Re^i2 - Re A3 + Q2) + ^3^12 ^ 



4Re^23 + 2 JA 



12 



12 



*We drop here the label "peak" for U"^, A''^ , and 5, for simplicity. 



(51) 
(52) 



The denominator of eq. (52) becomes under some conditions that are related to the value 
of the MNS matrix elements, 

2Re^23 + ■^^12 = , (53) 

where .4.23, J, and Aj2 denote some fixed value of them. By using the mixing angles 9ij 
and CP phase 0, written as 

2 

A12 = (cos (f) — tan 9i2 tan 6*23 sin 6*13) . (54) 

sin0 ^ ' 

Since 5 does not diverge at the first peak, the numerator also has to be 

(Re^l2 - ^-6^23 + "2) + «3^12 = , (55) 

under the same condition. Using eq. (53), 02 and Og become 



-YleAi2-, a^ = —-J. (56) 



Finally, we obtain 



2 fRe^,2 - Re^,2 + ReXg) - JA12 



and 

/i(Epeak) = 4 \U^,\^ \U,f - 4Asin(/. A12 - ReA2AL + 0(A?2) • (5^ 

Because the order of A^j is less than 1, the shift of the first peak for the transition prob- 
ability is not large. Equation (58) shows that the first-peak of the transition probability 
with ^MNS = 90° is smaller than that with (5mns = 270°. Since the value of ReAi2 is 
negative for p ~ 3.0 (g/cm^), the first-peak of transition probability for the CP-conserved 
case is slightly larger than that with A12 = 0. These features remain unchanged for the 
inverted hierarchy. But the shifting direction of the first-peak is different between the 
normal hierarchy and inverted one. 



5 summary 

In this article, we have estimated the three-fiavor effect for the determination of the mass- 
squared difference and the mixing angle from the survival and transition probability of 
u^. Prom both probabilities, the larger mass-squared difference is changed by the three- 
flavor effect. The order of magnitude of the difference between the larger mass-squared 
difference Sml^ and the mass-squared difference of the two-flavor analysis 5dip,peak is not 
only proportional to (5m^2 but also to the MNS matrix elements. We also flnd the CP 



phase dependence for the transition probabihty. If there is no three-flavor effect for i/^ 
survival and transition probabilities, the first-dip energy of z/^ survival probability and the 
first-peak energy of z/^ transition one are as same as that from A13. However, each of them 
is different from the value of A13, which are shown from eq. (30) and eq. (57), because 
of three-flavor effect. This means that the value of ^m^g from z/^ survival probability is 
slightly different from that from z/^ z/g transition one. The order of them are not so 
smaller than that of the expected error for the Smf^ in the future long base-line neutrino 
oscillation experiment. These results are useful to estimate the three-flavor effects for the 
value of the parameters which are obtained from the numerical analysis. 
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